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It has been recently shown that it is possible to excite gravitinos in an expanding background 
due to time-varying scalar field oscillating at the bottom of the inflationary potential. The two 
components of the gravitino, namely helicity 1/2 and helicity 3/2 are excited differently due to 
the presence of different time-varying mass scales in the problem. In this paper we analyse the 
production of both the helicities in multi-chiral scenario, in particular focusing on a general model 
of hybrid inflation. Fermion production in hybrid models is very much different from that of the 
chaotic models discussed so far in the literature. In this paper we give a full account of gravitino 
production analytically and numerically. It is noticed that the creation of gravitinos does not take 
place in the first few oscillations of the inflaton field, rather the production is a gradual and delayed 
process. It takes roughly 30 — 40 oscillations to build up the production and for the saturation to 
take place it can even take longer time, depending on the model parameters. We give an estimation 
of the reheat temperature and a brief discussion upon back-reaction on the gravitino production 
which could change its abundance. 



I. INTRODUCTION 

Low energy effective N = 1 supergravity is a predictive 
theory which could predict an inflationary potential 
flat enough to provide adequate density perturbations Q . 
So far, such viable inflationary models were constrained 
from observations by fixing the height of the potential, 
which essentially determines the amplitude of the COBE 
normalization. The first and the second derivative of the 
potential determines the tilt in the power spectrum, and 
the Yukawa couplings of the inflaton to other particles 
determines the reheat temperature of the Universe. The 
higher the coupling constant is, the higher is the tem- 
perature of the thermal bath and so the creation of the 
gravitinos from the collisions or the decay of other par- 
ticles. The gravitinos decay very late and depending on 
their mass their life time could be long enough to disrupt 
the synthesis of light elements, via hadronic shower or 
by altering the entropy density of the baryons. Due to 
these reasons there is a strong constraint on the reheat 
temperature, for a review see Ref. ||. 

However, there is also a non-thermal phase of the Uni- 
verse, just after the end of slow-roll inflation when the 
scalar field begins to oscillate coherently at the bottom 
of the potential. During this era an explosive produc- 
tion of particles, both bosons Q and fermions (|. |] , may 
take place due to the non-perturbative decay of the in- 
flaton to other fields. It has also been shown that it 
is possible to create super-massive bosons and fermions. 
However, fermionic production is always saturated by the 
Pauli blocking. As a matter of fact creation of heavy 
non-thermal bosons can be a good candidate for weakly 



interacting massive particles, known as WIMPS ||. The 
decay of super-massive bosons can explain the ultra-high 
energy cosmic rays 0. Super- massive fermions can be 
used in lcptogenesis, mainly from the decay of right- 
handed neutrinos into Higgs and leptons, which explicitly 
violates the CP conserving phase (8| . It is worth mention- 
ing that the non-perturbative technique of decaying infla- 
ton to other particles has given a new paradigm shift in 
understanding the hot big bang universe from the ultra- 
cold inflationary regime. 

Recently preheating in the context ofglobal supersym- 
metric theories has been considered ||,|l0| and, as a nat- 
ural extension it has been necessary to consider a lo- 
cal version of the supersymmetric theory and discuss the 
non-perturbative aspects of particle production and their 
consequences to nucleosynthesis. The local version of su- 
persymmetry, known as supergravity, naturally accom- 
modates the graviton and its superpartner the gravitino, 
a spin 3/2 particle. Quantization of spin 3/2 particles 
in the presence of an external background is plagued by 
consistency problems. It has been known for a long time 
that quantization of spin 3/2 particles in scalar, electro- 
magnetic, or gravitational backgrounds can give rise to 
acausal behaviour pdj |. Supergravity is the only set-up 
where such problems do not occur, provided the back- 
ground fields also satisfy the corresponding equations of 
motion jlj] . Nevertheless, the complicated form of the 
Rarita-Schwinger equation makes it extremely difficult 
to extract any explicit results even in a simple back- 
ground . The problem was first addressed in Ref. Q , 
where the authors quantized spin 3/2 particles in a non- 
vanishing cosmological background, almost two decades 
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ago. 

The slightest generalization of quantizing spin 3/2 has 
been done very recently in literature, in Ref. Jl5[ ] . The 
authors have extended the calculation of quantizing spin 
3/2 in presence of a time- varying homogeneously oscil- 
lating scalar field in a cosmological background. This is 
the first result where the non-perturbative decay of the 
inflaton to gravitinos during preheating has been taken 
into account. The authors have explicitly shown the pro- 
duction of a particular helicity, the ±3/2 components 
of gravitino, in a particular new-inflationary type model 
]16j . It has been noticed that the non-perturbative result 
can give rise to a larger abundance compared to the per- 
turbative decay. The gravitino to photon number density 
has been found to be n^^/n^ ~ 10~ 12 Jl^]. This abun- 
dance is 3 orders of magnitude larger than the thermal 
abundance for a reheat temperature 10 5 GeV j|. Such 
over production of gravitinos has been the first proof of 
the non-thermal production of gravitinos with helicity 
±3/2, which demands constraining the reheat tempera- 
ture in any supergravity motivated inflationary model. 

However, a massive gravitino has 4 degrees of free- 
dom, the other two degrees of freedom due to the helic- 
ity ±1/2 components of the gravitino. The production of 
helicity ±1/2 gravitinos is directly related to the problem 
of super- Higgs mechanism in supergravity models, which 
was studied in Ref. [ fl4| in the context of a non- vanishing 
cosmological constant, and in Ref. 17 with a vanishing 
cosmological constant. In the presence of a time- varying 
scalar field there is an additional source of supersymme- 
try breaking via the non- vanishing time derivative of the 
homogeneous scalar field. This plays an important role 
in the context of cosmology when the scalar field is rec- 
ognized as the inflaton, oscillating coherently at the bot- 
tom of the scalar potential. Due to the presence of such 
a field, supersymmetry is always broken at the minima of 
the potential and the initially massless gravitino which 
possesses only the helicity ±3/2 components "eats" the 
goldstino to gain the other ±1/2 components. At this 
point, one may wonder how to generalize the super-Higgs 
mechanism in such a scenario. In fact the problem turns 
out to be quite complicated and it has been addressed in 
two seminal papers Jl^Jl^]. These papers also study for 
the first time the production mechanism of the helicity 
±1/2 components of gravitino ( see other papers in the 
similar context |20| ) . In all these papers the authors have 
assumed the existence of a unitary gauge, where the phys- 
ical Lagrangian is free from the goldstino field. We will 
explicitly show that it is possible to choose such a gauge, 
where the gravitino equation of motion is free from the 
goldstino field. Our calculation is valid for more than one 
chiral field as well. Though we shall give a proof for the 
-F-type supersymmetry breaking, this can be extended to 
D-type supersymmetry breaking also. 

It has been noticed in Refs. Jl8|,[l9| that the production 
of the two helicity states are completely different. The 
helicity ±1/2 components are produced copiously com- 
pared to the helicity ±3/2. In the case of helicity ±3/2, 



conformal invariance is broken due to the presence of 
a time- varying gravitino mass, which is usually Planck 
mass suppressed, whereas for helicity ±1/2 the breaking 
of conformal invariance is related to the presence of a 
massive Goldstone fermion [ ^8[Jl9| ], whose time- varying 
mass is not suppressed by the Planck mass. Moreover in 
the high momentum limit the helicity ±1/2 gravitino be- 
haves like a fermion, the goldstino, as it is stated by the 
equivalence principle |2^j2^]. This has been studied in a 
single-chiral field scenario, where the source of conformal 
breaking can be directly related to the mass of Goldstino. 
Unfortunately in the multi-chiral field scenario the quan- 
tization scheme becomes more involved, and the relation 
between conformal breaking and Goldstino mass is not so 
straightforward. The situation has been briefly discussed 
in Ref. pl[ | , and an attempt of a perturbative scheme has 
been suggested. However, it would be nice to discuss a 
non-perturbative scheme. Our paper fills that gap and, 
as we shall see, we can discuss non-perturbative produc- 
tion of the helicity ±1/2 gravitinos in the multi-chiral 
case. 

The best example to study the multi-chiral field sce- 
nario is in the context of a general class of supersymmet- 
ric hybrid inflation model |^,[n|. There are essentially 
two scalar fields, one is responsible for inflation, and the 
other field is responsible for the phase transition which re- 
sults in terminating the inflationary era. Unlike the non- 
supersymmetric version of hybrid inflation model |p3| , 
the supersymmetric version considered in this paper has 
only one coupling constant in the potential. This leads 
to a single natural frequency of oscillations. This gives us 
an ample opportunity to use techniques to explore grav- 
itino production in a similar spirit as in the case of a 
single-chiral field. Fcrmionic production in this case is 
very much different compared to that of a quadratic in- 
flationary potential. In the hybrid models the effective 
mass term for the fermions is always positive and as a 
result the production can never be completed in just a 
few oscillations, rather the occupation number gradually 
increases and depending on the model parameters the 
rate of production could be slowed down significantly. 
For example, in order to reach the Fermi saturation it 
may be necessary more than 100 oscillations. Due to 
such a slow rate of production, the issue of back-reaction 
becomes important. It is very likely that simultaneous 
non-perturbative production of bosons can change the 
picture quite significantly. In some sense, hybrid model 
can be considered to be the safest of all the supergravity 
motivated inflationary models, because the gravitino pro- 
duction can stop due to back-reaction effects coming from 
other newly created bosons or fermions. In other models, 
such as in the chaotic models, where most of the particle 
creation takes place in the very first few oscillations, the 
issue of back-reaction hardly plays any significant role. 

The layout of the paper is as follows; in section 2, 
we establish the supergravity Lagrangian and discuss the 
gauge fixing mechanism for the multi-chiral field scenario; 
this is a mere generalization of the gauge usually used 
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in pure gauge theories as well as an alternative way of ex- 
plaining the existence of a unitary gauge in supergravity 
theories [pi 25 2^]. In section 3, we discuss the quantiza- 
tion procedure for the gravitinos. In section 4, we briefly 
describe a general class of supersymmetric hybrid infla- 
tionary models we use. Analytical results and discussions 
based on our numerical results are presented in section 
5. The implications of these results on the reheat tem- 
perature are discussed in section 6. We give a detailed 
discussion on fermionic creation in hybrid model in the 
appendix. 



II. GENERAL SUPERGRAVITY LAGRANGIAN 
AND GAUGE-FIXING 

In this section we describe the supergravity La- 
grangian. For the sake of brevity, and for our purpose, 
we concentrate upon the chiral-supermultiplet which con- 
tains the bosonic part, fermionic part and the interac- 
tion terms between the fermions. We consider minimal 
Kahler potential with K = 4> l 4>i, where the scalar fields 
(j>i are taken to be real. The superpotential is denoted 
by W, and the Kahler function is given as usual by 
G = K + ln|M^| 2 . The choice of minimal Kahler po- 
tential also ensures G* = <5* . The total Lagrangian is as 
follows 

e^suGRA = I + dubVp + e G (3 - G.G 1 ) 



i g G/2 (Qij + Qi G j\ 



XiXj 



(i) 



where indices represent species of the chiral 

multiplets and qt = j x d\. The derivatives — 
— ^Lo^mn& mn are covariantized with respect to grav- 
ity, where w pm „ is the standard spin connection [p6| . 
The gravitino field is represented by Xi represents 
the fermions and e is the determinant of the vier- 
bein e£. We have also used a» v = |[7 M ,7 !/ ]. The 
complete Lagrangian would also contain the fermion 
Yukawa couplings, four fermion terms, and numerous 
non-renormalizable terms, details can be found in the 
literature M. Here we have neglected the torsion term 
in the covariantized derivative, assuming that the grav- 
itino production is small, so their back reaction can be 
neglected. However, we mention that consistency of su- 
pergravity theory would demand the presence of the tor- 
sion term in the covariantized derivative. These terms 
are suppressed by the Planck mass square and we can 
safely neglect them in the linearized gravitino equation 
of motion. The Planck mass Mp has been taken to be 
unity and the convention shall be retained, unless other- 
wise stated. 



The above Lagrangian is invariant under local super- 
symmetric transformation laws |26f| . For spontaneous su- 
persymmetry breaking to occur at least one of the field 
vacuum expectation values should be non-zero. In par- 
ticular, for _F-term type breaking of local supersymmetry 
one requires 

-i - e G/2 GiO ± , 



<ofe|o> 



(2) 



where £ is the infinitesimal Grassmann-odd parame- 
ter. The right-hand side of Eq. (|^) has two explicit 
terms which can break local supersymmetry. The sec- 
ond term is the usual F-term of the scalar field whose 
non-vanishing vacuum expectation value induces super- 
symmetry breaking. The first term will give a non-zero 
contribution in case of a time- varying scalar background 
field. This will be the situation if we identify the back- 
ground fields 4>i with the oscillating inflaton fields. In 
this case local supersymmetry is always broken during 
the oscillations. 

The goldstino can be identified as usual from Eq. (||) 

V = 0*X l , (3) 

where, 

9t = i ficfri ~ e G ' 2 Gi , (4) 

and we follow the notations introduced in Ref. [ pi) . We 
are only interested in homogeneous scalar fields (f>i which 
are solely function of time. The goldstino is then "eaten" 
by the gravitino in local supersymmetric theories. This 
also ensures that the gravitino gains the helicity ±1/2 
components other than the ±3/2 components, and be- 
comes massive. In the high energy limit it is possible to 
relate the helicity ±1/2 components of the gravitino to 
the goldstino via the equivalence principle [ p7| , plp^ ] . In 
the limit when Mp — > oo, the helicity ±1/2 components 
retain the memory of the goldstino contribution and this 
is the reason why the two helicities behave differently and 
have different production rates Jjl|[l!| . 

As it can be realized by inspecting Eq. ([!]), the grav- 
itino is coupled to the fermions through the following 
mixing terms 



7Y- 



(5) 



However, it is possible to get rid of these mixing terms by 
adding a gauge-fixing term to the Lagrangian as it has 
already been discussed in various places p^^5| , p2[ . In 
this section we mainly concentrate on the mere existence 
of such a gauge- fixing term in supersymmetric theories. 
Regarding this, we extend the previous calculation made 
on gauge-fixing for a single-chiral fielcj^ |22] . 



*In Ref. |22|, the authors have only considered the equation 
of motion for the goldstino, because they were more interested 
in establishing the high-energy equivalence between helicity 
±1/2 gravitinos and goldstinos. 
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Before we move onto specifying the gauge, let us in- 
troduce the projection operators, which we need to use 
latex on MM. 
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'3 ! 



V 1 — — 6 

The modulus of 9i is given by 

0*6 = e G G i G l + ^'(t)0i(t) = p + 3e G : 



(6) 
(7) 

(8) 



where derivative with respect to time is denoted by dot, 
p — <p l 4>i + V is the energy density, and the scalar poten- 
tial is V = e G (G l Gi — 3). Accordingly, the fermion Xi 
can be split into two components by using the projection 
operators Eqs. (||) and (j7|) 

Xl = x * + We v ' ^ 

Xt = T>& ■ (10) 

Now, with the help of Eqs. (H|) and (|J) the mixing terms 
between the gravitino and the chiral fermions can be re- 
cast as 



V2 



where we have used 



(ii) 



(12) 



In a similar way, one can also reduce the complex con- 
jugate part of the mixing terms. Eq. (|ll]) tells us about 
the direct coupling of the gravitino to the goldstino. The 
equation of motion for the gravitino ip^ that follows from 
the Lagrangian Eq. ([!]) is then 



fin 



V2 



-?7 = 0. 



(13) 



with an explicit term depending on the goldstino field due 
to the mixing. However, we notice that the contribution 
of the mixing term in Eq. ( pi] ) can be canceled by adding 
to the Lagrangian the following gauge-fixing term 



i(F pF , 



(14) 



where p = j x D\, F = F'^oi and the gauge-fixing func- 
tion is given by 



F(i>, V ) = 



We 



1Pu 



1 1 



(15) 



This is a mere generalization of the gauge used in pure 
gauge theories. Similar gauge-fixing has been initially 



introduced in Ref. p9] in the static case, where the scalar 
field has been taken to its value at the minimum of the 
potential. Therefore, once we fix the gauge, the equation 
of motion for the gravitino is completely free from the 
goldstino. The limit £ — > corresponds to the unitary 
gauge, and it implies — > 0. This is equivalent to demand 
that no goldstino be present in the physical spectrum. 
From here onwards we will work in the unitary gauge. 

Using the gauge-fixing condition, F(i/j,r)) — 0, and de- 
manding £ — > 0, the final equation of motion for the 
gravitino, namely Eq. (13) in the unitary gauge can be 
written as 



,-1 Aivpa 



i e G /V>7l^ = 0. (16) 



Here the mass term for the gravitino e G / 2 depends on 
time, due to presence of oscillating background scalar 
field <pi(t), whose dynamics we shall discuss later on. 
From now onwards we express the mass term as m(t) 



m(t) = e G ' 2 = e K ' 2 \W\ 



(17) 



where W is the model dependent superpotential. De- 
tailed discussion will be given in the coming sections. 
The above demonstration of removing the goldstino de- 
pendence from the gravitino equation of motion in the 
unitary gauge suggets that it is possible to generalize 
gauge for a multi-chiral time-varying scalar background. 
Imposing the unitary gauge simplifies the gravitino field 
equation of motion in general and now we will be in- 
terested in quantizing gravitino field in a cosmological 
background. 



III. DIFFERENCE BETWEEN HELICITY 1/2 

AND 3/2 

We have seen in the preceeding section that we can 
recognize the goldstino component which couples to the 
gravitino field in a scenario when the dynamics of the 
background scalar field is also taken into consideration. 
By using the gauge-fixing term in the Lagrangian, we 
have shown that it is possible to cancel the gravitino- 
goldstino coupling term appropriately. The final equa- 
tion of motion for the gravitino field, thus free from the 
goldstino, can be used to study the gravitino production 
in a dynamical background dominated by an oscillating 
scalar field. In this section we do not attempt to re- 
derive the equations of motion, which have already been 
discussed in several papers 15 l^,|l^,^l[, rather we dis- 
cuss few issues and the main equations. By studying 
the equation of motion for the Rarita-Schwinger field, 
one notices that there is a free index left, which in prin- 
ciple can be contracted by at best two possible ways, 
say 7^ or Z? M , thus giving rise to two constraint equa- 
tions for the whole system. In presence of a cosmological 
constant, the equations of motion for both the helicitics 
look alike, with two simple constraint equations, namely 
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7 M Vv = 0, D^ijj^ = 0. However, this is not correct in any 
arbitrary gravitational background. These constraints do 
not hold true for the helicity ±1/2 case in an oscillating 
scalar background [|l8Ul9| , even though these constraints 
continue to hold for the helicity ±3/2 case in the same 
oscillating background jL5|. The helicity ±1/2 compo- 
nents gain an effective mass during the oscillations of the 
inflaton, but the helicity ±3/2 do not seem to see the 
effect of curvature at all. This gain in mass is purely due 
to presence of a non-trivial background curvature, which 
suggests that the different helicities couple to gravity dif- 
ferently. 

We follow the notations of Ref. in order to write 
down the equations of motions for helicity 1/2 and 3/2 
for a momentum mode (0, 0, 0, |fc|) projected along the z 
direction. 

(i'fdo + »~7° - ma + k^ i/> 3/2 (t, a:) = , (18) 
^7 ^,+t~7 -mo + *G7 3 ^ 1/2 (r J a:) = 0. (19) 



The equations have been written in conformal time dr = 
dt/a, where a is the scale factor. Prime denotes derivative 
with respect to r, and m is the gravitino mass Eq. (17). 
Here V'1/2 and V'3/2 ar e Majorana spinors which would 
correspond in the flat limit case to the 1/2 and 3/2 he- 
licity states respectively (for details see Ref. 119]). The 
matrix G in Eq. (|l9| ) can be expressed in terms of A and 
B functions |J,[ij 



G = A + ij°B 



p — 3m 
p + 3m 2 



■ «7 



2m' a- 1 
p + 3m 2 



where p and p are denoted by 



p = £|4| 2 + 7(^), 

i 



with 



V{fa) = e K (\diW + faW\ 2 - 3|Wf 



(20) 

(21) 

(22) 
(23) 

(24) 



Here again dot means time derivative with respect to 
physical time. The potential energy Eq. (p4| ) is given in 
terms of the Kahler potential K and the superpotential 
W. In the limit when fa -C 1 (in units of the Planck 
mass), A and B can be expressed in a simpler form 



G = A + i 1 °B = P - + i 1 *™. 

P P 



(25) 



It is important to point out that in general \G\ 2 = A 2 +B 2 
is time dependent, and \G\ 7^ 1. Only in the case of 



a single-chiral field \G\ = 1. For multi-chiral field sce- 
nario, in particular in the case of a supersymmetric hy- 
brid model, \G\ departs from 1. This makes the quantiza- 
tion scheme slightly more involved than simple scenarios 
where \G\ = 1. To proceed with the quantization we 
redefine G in terms of the conformal time 



G = A{t) + z 7 °S(r) - el adr e 2ll ° I ^ 



(26) 



where the coefficient in front of the overall phase rep- 
resents \G\. We concentrate upon the helicity 1/2 case, 
helicity 3/2 being a simpler generalization of that. We 
expand V'1/2 m terms of the mode functions 



1/2 



-5/2 



d 3 k 
(27T) 3 / 2 ' 



-ik.x , 



-ry J fidr 



(u r (r,k)al +v r (T,k)b[ 



k) ' 



(27) 



where, v r (r,k) = u r (r, — k). The spinor u r (r,k) satis- 
fies the following equations of motion 



z fim f[U± + ik\G\u 



m e g = ma + p , 



(28) 



where u T = (u+,u_), and for the gamma matrices we 
have used the representation in which 
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1 
-1 



.7 



1 

-1 



(29) 



It is possible to write down a second order differential 
equation from the set of equations in Eq. (Eq) . 
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u ± 



&4 + ('' J|r;|J 



m cS ± i 



H 

\G\ 



Weff 



Eq. (£0D 

e -J <*/2dT 



T im' cS ) u± = . (30) 
can be further reduced by redefining u± — » 



a 1 



(k 2 \G\ 2 + n 2 Titt')u ± = 0, 
il = m c ff ± t — 



or, by redefining a new time in Eq 



(31) 



d 2 u± 
dz 2 



+ k 2 + 



TOeff 

\G\ 




where m c g is defined in Eq. ( |28| ) , and the equation is for- 
mally analogous to the evolution equation for a spin-1/2 
fermion in a time- varying background. It is important to 
notice that all the three equations Eq. ( |30| ) , Eq. (|3l]) and 
Eq. ([52]) are equivalent, expressed in different forms. For 
our numerical results we have used Eq. (30), and for our 
analytical treatment we can consider any of these three 
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equations. In the next section it will become clear that 
\G\ < 1 for any general model, so that the factor \G\ 2 
in front of the momentum squared cannot give rise to 
any acausal behavior. Nevertheless, we note that par- 
ticle creation could take place due to time variation in 
\G\. However, from Eq. ( |32| ) it is clear that this requires 
a non- vanishing effective mass term m c ff, even if it is 
merely a constant. Particle production takes place due 
to the breaking of the conformal invariance. In our case, 
the violation of conformal invariance is not only due to 
presence of the gravitino mass but also due to /i, which 
can be related to the goldstino mass (l^Jl^,^l| . 

In order to evaluate the occupation number we first 
evaluate the hamiltonian 

H(t)= I d 3 kJ2 E k(r)(ata r - b r bt) + F k (T)b r a r 

J r 

+F* k {r)albl), (33) 

in which the momentum k is along the third axis, and 

E k = 2k\G\Re(u* + u_) + m cS {\u_\ 2 - \u+\ 2 ) , 
F k = 2m e ffU_)_M_ + fc|G|(u^_ — u_) , 



El 



IF, 



2 1.2 



'off 



(34) 



The Hamiltonian can be diagonalized with the help of a 
Bogolyubov transformation. The new set of creation and 
annihilation operators which diagonalize the Hamiltonian 
can be defined as 



a(fc,r) =a fc (r)a(fc) + /3 k { T )b\-k) 
b\T) = -f3i(T)a(k)+at(T)tf(-k) 



(35) 



where a k , and (3 k are the normalized Bogolyubov coeffi- 
cients, 



Ok 

(3 



Ek 



■ to 



F* 



2u(uj + E k ) 



QJ - E k 
2uo 



t off 



\G\ 2 k 2 



(36) 



Now, the time dependent occupation number can be writ- 
ten in terms of the vacuum expectation value of the num- 
ber operator 



n(r) = (0|7V|0) = 



1 



7r 2 a 3 (r) 



dkk 2 \(3 k \ 



(37) 



In order to solve Eq. (31j), one needs to specify the bound- 
ary conditions. Usually they are defined such that at the 
beginning, when t — > 0, we have \Pk\ 2 — 0, and the oc- 
cupation number n(0) = 0, suggesting that there is no 
particle density at the initial time: 



«±(0) = 



u' ± (0) = Tim e gu ± {0) + ik\G(0)\u T (0) . 



(38) 



Now, we have all the tools necessary to study the produc- 
tion of gravitinos in the hybrid inflationary model. Next 
we describe the hybrid model and how to estimate the 
occupation number. 



IV. HYBRID INFLATION 

The hybrid inflation potential can be derived from the 
following superpotential 



W = M{N 2 - N 2 ) , 



(39) 



where <f> plays the role of inflaton. During inflation the 
other field N is trapped in its false vacuum, N = 0, 
while the <f> field rolls down towards the critical value 
determined by the value of N at the global minimum, 
A^o, and the coupling constant A between <f) and N. At 
this point the N field rolls down from its zero value to- 
wards the global minima and begins to oscillate around 
Nq, while (j> oscillates around zero. This will enable the 
preheating phase of the Universe. During this period an 
effective potential for the fields 4> and N can be derived 



V= \W^\ 2 + \W N \ 2 , 

= \ 2 {N 2 - N%) 2 + 4\ 2 (j) 2 N' 



(40) 



where the subscripts denote the derivative of W with re- 
spect to the fields. The superpotential in Eq. ( |39| ) also 
ensures a non-vanishing constant vacuum energy during 
inflation V(0) — \ 2 Nq. We should also mention that A 
and Nq act as free parameters of the model, but they are 
constrained to some extent from the COBE normaliza- 
tion and the tilt in the power spectrum |L0| 



AA^o re 1.27 X 10 15 |?7*|GeU, 



(41) 



where we have taken[j] |?/*| sa 0.01 in our analysis, which 
is a reasonable assumption in order not to generate a 
sharp tilt in the power spectrum of the density pertur- 
bation during inflation. The present constraint on the 
spectral index is \n — 1| < 0.2 0. Hybrid inflationary 
model derived from such a superpotential is known as F- 
term hybrid inflation. Slightly different version of hybrid 
inflation popularly known as D-term inflation pp| can be 
derived from the Fayet-Illiopoulus term appearing from 
an anomalous U{1) symmetry, which could provide the 
necessary potential energy during inflation. Whatsoever 
be the cause of such a potential, our argument for the 



t?7* is one of the two slow-roll parameters, usually defined 
in inflationary cosmology as ry» = (mj>/(87r))(V"/V(0)); the 
"star" index means that it is evaluated at least 60 e-foldings 
before the end of inflation. This should not be confused with 
the goldstino; we have already defined goldstino with the same 
notation earlier. 



G 



gravitino production is quite generic and will not depend 
upon a particular origin of the vacuum energy. 

Due to presence of a single mass scale XNq in the 
model, which is related to the supersymmetric breaking 
scale during the inflationary era, we will have a single fre- 
quency of oscillations during the preheating phase, and 
effectively a single scalar field oscillating. In our case this 
can be taken as N, related to the other field <fi by 



N - N(t) 
V2 



(42) 



By solving the equation of motion for N(t) |^|TC]] we get 
the approximate solution when the field is around the 
bottom of the potential 



N(t) 
Nn 



1 + S(t) cos(iV) , 



(43) 



where v$ — 2XNo provides the natural frequency of the 
oscillations, and £(t) is the amplitude of the oscillations 
decreasing in time as ~ 1/t [. The above expression 
is valid when \N{t)/N - 1| < 1/3, with E(0) ~ 1/3. 
What makes the hybrid scenario interesting is that N(t) 
never vanishes and that causes the fermionic creation to 
be completely different than the usual chaotic type infla- 
tionary potentials. 

Now, with the present knowledge we can evaluate |G| 
with the help of Eq. @ 



|G| 



a { P -2Vf , 4|PU| 2 



P z 



'N<, 



4A 2 iV 2 • „, N9 

= 1 j^-\N\ (N — N ) , 

P 



(44) 
(45) 
(46) 



where the subscript in Wa, , Wn means derivative with re- 
spect to the corresponding field, and dot denotes as usual 
physical time derivative. The last equation has been writ- 
ten using Eq. @, knowing that \cf>\ 2 + \N\ 2 = 3/2|iV| 2 . 
Here it is important to notice that the departure from 1 
is quite obvious, even though there is effectively a single 
scalar field oscillating. In fact, we can easily general- 
ize Eq. (E^) to any arbitrary number of scalar fields fa 
contributing to the energy density. Using Eqs. (pfj|-p4[), 
without assuming fa <C 1, we get 



|G| 2 = 1- 



K 



(p + 3m 2 )' 



^ V,\Vo, VjWo, 



(47) 



i<j 



where T>iW = Wi + K{W, and the subscript i denotes 
derivative with respect to the homogeneous scalar field 
fa. This shows explicitly that |G| < 1. 



'•The actual decrease in amplitude over an oscillation period 
depends on the ratio H/v$, H being the Hubble parameter. 



V. ANALYTICAL ESTIMATION OF 
OCCUPATION NUMBER 

A. Helicity 1/2 

It is possible to analytically estimate the occupation 
number, and the number density n(t) for both the helic- 
ities. Assuming that |/3fe| 2 ~ 1, for a given momentum 
k, our task reduces to estimate the cut-off momentum 
k. For this we need to solve Eq. ([jl|), which will then 
lead us to study the occupation number for the helicity 
1/2 gravitinos in hybrid model. In the case of helicity 
1/2 states, the contribution of |G| 2 can be important, 
and we need to evaluate it before we could estimate the 
abundance of helicity 1/2 states. Later on we will discuss 
the abundance of helicity 3/2 states, where |G| 2 = 1. To 
carry out our calculation we need to know the dominant 
contribution to appearing in Eq. (pll), 



where 



O = ma + p + i — , 



j¥_ A \G\' 
~2B + 2B~\G\ 
\GY 



\G\ • 



(48) 

(49) 
(50) 



With the help of Eq. Q and Eq. (g) it is possible to 
evaluate all the terms in Eq. fliq). Here we simply quote 
the final results in terms of the physical time: 



m{t) 

|G| 2 

\G\ 
\G\ 
A 

2B ~ 2V2 
A \G\ 



4 V / 2A 2 M 2 



(l + cos(2^i))+C(l/i 3 ), 



1- - sin 2 (2^) + 0(l/t), 



(51) 
(52) 



va, , , / sin 2 (2zAj,i) \ 
-^sin(4^)(l + \-±l) +0 (l/t), (53) 



3^ ( 1+ 5m cosM \ +0{1/t 2 ) 



(54) 



2 B M = ^ cosW)(l + 



sin 2 (2r/0t) 



-0(l/t). (55) 



It is apparent that the only term which dominates £l(t) in 
Eq. ( pq ) is due to Eq. (p4[). The mass term m(t) is sub- 
dominant due to the Planck suppression ( here we have 
explicitly written the Planck mass). The amplitude of 
the oscillations from Eq. ( |55| ) is one order of magnitude 
smaller compared to that in Eq. (|54|). As a first approx- 
imation we can consider A/(2B) to be the effective mass 
term for the helicity 1/2 states. 

A general discussion of production of massive fermions 
in an oscillatory background is given in the Appendix. 
We will use the results quoted there in order to esti- 
mate the occupation number for the helicity 1/2 grav- 
itinos. By taking Q(t) » -A/(2B), and |G| 2 « 1, we 
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notice that Eq. ([H]) mimics Eq. ( |Al[ ) in the appendix 
for a time-varying mass which is denoted by fl(t) in our 
present situation. First of all we notice that the effec- 
tive mass term A/(2B) never vanishes at any point of 
the scalar field oscillations. With £(f) ~ 1/3, the effec- 
tive mass is always positive. This situation is ty pica l of 
any hybrid model. Comparing Eq. ( |54| ) with Eq. dA2| ) in 
the Appendix, we have i>$ — m^, mx — 3i/^/(2v2) and 
g(p(Q) = 5i/^/(8v2). It is evident that mx > 50(0), and 
as a result m(t) is always positive, quite similar to the 
case of bosonic production, where the mass term appears 
always as squared. However, when studying fermionic 
creation for a chaotic type potential there is a possibility 
of having m(t) vanishing at some point during the oscil- 
lations B . In this case fermion production takes place in 
the first few oscillations, since the adiabatic condition is 
violated maximally when the inflaton field passes through 
the point where the effective mass vanishes, and it is pos- 
sible to create very heavy fermions in the process. In hy- 
brid models and for the helicity 1/2 gravitinos, because 
the effective mass remains always positive, the adiabatic 
condition is broken when the effective mass reaches its 
minimal (non zero) value, and the mass of the fermions 
created never exceeds the mass of inflaton. In this case 
it is also important to note that the amplitude of the 
oscillations decays slowly and it takes roughly 20 — 30 
oscillations to make any significant change in the initial 
amplitude. As a result the production of fermions takes 
place at each and every oscillation, and since the degree of 
violation of the adiabaticity is much weaker compared to 
the former, the production process takes a longer time to 
saturate the Fermi band. This is quite evident from our 
numerical results, see Fig. (fil) and Fig. Q). The choice of 
model parameters in Fig. (ffl) leads to a higher inflation- 
ary scale compared to that of Fig. (||). However, the rate 
of production is exactly the same in both the cases. The 
reason is that the occupation number essentially depends 
on Vfj,, which is exactly the same in both the models we 
have considered. This can be vividly seen by comparing 
the number of peaks in the occupation number. 

In Fig. U it is easy to see the particle creation taking 
place in bursts and in regular intervals. The peaks re- 
mind us the violation of the adiabatic condition. This is a 
perfect example of particle creation in a broad resonance 
regime. In the narrow resonance, particles are always 
created throughout the evolution of the scalar field [Q, 
they are not produced in bursts like in Fig. |^. It is visible 
that the particle number remains constant for some time 
and then there is a sudden jump in the occupation num- 
ber. Physically one can understand the situation from 
the comparison with an interacting quantum field the- 
ory. Usually in field theory, at extreme past/future we 
consider the plane wave solution of the free equation of 
motion. However, in between extreme past and future 
the interaction is switched on. Here also as we can see 
that the occupation number remains constant for a while, 
which represents the ingoing wave and then the interac- 
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FIG. 1. The evolution of the classical field N(t), and the 
number density of helicity 1/2 gravitinos (scaled by a factor of 
10) has been depicted for the choice of parameters A = 10 -3 
and No — 2 x 10 16 GeV. The time scale is given in terms of 
the approximate number of oscillations of the fields in physical 
time, JVos = v^t/2n. The number density is given in units of 

tion switches on, which is depicted by a jump in the occu- 
pation number. The outgoing wave in this picture is the 
newly filled occupation number. In fact the Bugolyubov 
coefficients can be calculated by estimating the reflection 
and the transmission coefficients of the plane wave pass- 
ing through a barrier. Most of the important information 
concerning the production of helicity 1 /2 gravitinos can 
be extracted from Fig. ||. The production of helicity 1/2 
gravitinos builds up in each and every oscillation. The 
production does not saturate the Fermi level in the first 
few oscillations, but it takes several oscillations to reach 
the Fermi level. This behaviour is in stark contrast to 
the fermionic creation in a quadratic potential || , where 
the Fermi level is reached in a few oscillations. It is also 
noticeable that there is no stochastic behaviour in the 
occupation number, and it is always increasing until it 
reaches the Fermi level. The main reason is that the 
effect of expansion is felt very slowly in hybrid models, 
when compared to the typical period of oscillation, that 
is, Hq /v,), w No /Mp <C 1 , where Hq is the value of the 
Hubble constant at the end of inflation. Therefore, the 
amplitude decreases very slowly in most of the cases, al- 
most adiabatically. There is also a slight change in fre- 
quency, as it can be noticed from the plot (for details see 
Ref. |10) ) . However this small change in frequency is not 
going to affect our analytical estimation. 

We have plotted the spectrum for the helicity 1/2 grav- 
itinos in Figs. (|3|) and for two different set of model 
parameters. In both the cases Fermi-level is saturated 
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FIG. 2. The evolution of the classical field N(t), and the 
number density of helicity 1/2 gravitinos (scaled by a factor 
of 10) has been depicted for the choice of parameters A = 1 
and No = 2 x 10 13 GeV. The number density is given in units 
of We note that the rate of production depends mainly 
on and as a result the production rate is similar to that 
of Fig. (§). 



for fc max ~ v&- This confirms our a nalytical study in 
the preheating section, see Eq. ( A14 ), which reduces to 
numerically observed fc max for rax ~ m<f>- The number 
density of helicity 1/2 can be obtained from Eq. (|37 



"1/2 



1 

47T 3 



d 3 ka(k) 



h 3 



37T 2 



3tt 2 ■ 



(56) 



Here fc max has been taken to be a comoving momentum. 
It is evident from Fig. (|l|) and Fig. (0) that the occupa- 
tion number grows gradually and saturates after many 
oscillations depending on the choice of A and Nq. 



B. Helicity 3/2 

We can follow similar arguments to evaluate the oc- 
cupation number for helicity 3/2 gravitinos. Notice that 
Eq. ( p^| ) for helicity 3/2 reduces exactly to Eq. ( |3l"| ) with 
\G\ = 1 and Q = m eS = ma, ||,p|,^9|Q. Hence, the 
single time-varying mass scale appearing in the prob- 
lem i s no w given by Eq. (^l|). Comparing Eq. (]5l|) with 
Eq. (A2), it becomes clear that the bare mass mx itself 
is time- varying for the helicity 3/2 case, but as we have 
noted before, the amplitude of the oscillations is almost 
constant, especially in the hybrid model we are inter- 
ested in. The parameters rax and g<p(0) are equal to 
each other, and both are Planck mass suppressed. But 
this does not mean that we cannot excite helicity 3/2 
gravitinos. The main point is that the effective mass 
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FIG. 3. The spectrum of helicity 1/2, when the number of 
oscillations is iV s = 60, for the model parameters A — 10~ 3 
and No — 2x 10 16 GeV. This choice of model parameters leads 
to -q, ~ 0.01. It is evident that the Fermi level is saturated 
for a cut-off momentum k ~ v$ ~ 2\Nq. 



always vanishes at some point during each and every os- 
cillation. As a matter of fact we can excite them pre- 
cisely due to this. What matters is the violation of adi- 
abaticity, and that takes place precisely at those points 
where m(t) vanishes. Production of helicity 3/2 mim- 
ics the first scenario we have discussed in the appendix, 
but now rax "C ra^ — Iv^. However, since the effec- 
tive mass ra(t) vanishes in each and every oscillation, so 
the gravitino production takes place in each and every 
oscillation. The spectrum of helicity 3/2 gravitinos is 
plotted in Fig. (g), for the model parameters given in 
Fig. (|3|). The spectrum preserves the essential features 
obtained before, but most importantly Fermi level never 
gets saturated and the production is extremely subdom- 
inant compared to helicity 1/2 case. Nevertheless, the 
helicity 3/2 gravitino abundance may also pose a strong 
bound on the model parameters and the reheating tem- 
perature, and thus it is necessary to study them as well 
]l5[ . Ou r task is to estimate fc max , and with the help of 
Eq. (|A|) we get 



,3/2 



2y/2 A 2 Ml 



(57) 



By taking M p w 10 18 GeV, and v$ w 10 14 GeV, we get 
fcmax ~ 0.11/0, which matches very well with our numeri- 
cal result. In Fig. (||), the spectrum peaks around O.3z/0. 
It is straightforward to estimate the occupation number 
for the helicity 3/2 gravitinos 



"3/2 



6V2^ 2 A 2 M 2 



(58) 
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FIG. 4. The spectrum of helicity 1/2, when the number of 
oscillations is iV os = 60, for the model parameters A = 1.0 and 
No = 2 x 10 13 GeV. This choice of model parameters leads to 
?7» ~ 0.01. It is evident that the Fermi level is saturated for a 
cut-off momentum k ~ v$ ~ 2AiVo . Important point to notice 
that this choice of model parameters leads to low inflationary 
scale and this gives a different spectrum than noted earlier in 
Fig. ©. 
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FIG. 5. Spectrum of helicity 3/2, when the number of 
oscillations is iV os = 60, for the model parameters A = 10~ 3 
and TVo = 2 x 10 16 GeV. This choice of model parameters leads 
to rj, ~ 0.01. It is evident that the Fermi level is not saturated 
in this case. Even though the production is small compared 
to helicity 1/2, it is certainly not negligible. Important thing 
to notice that the spectrum preserves its shape and peaks 
around 0.3f^. 



Here we have assumed (3k in Eq.(|37|) to be a constant, 
< |/3| 2 < 1, for the cut-off momentum k max . The actual 
value of P is difficult to estimate analytically, and we do 
not attempt to analyse it here. We can also estimate the 
abundance ratio for the two helicities 



"1/2 
"3/2 



2y/2\ 2 M£ 

WWW) 



(59) 



We should mention that in Eqs. (|57[-p9[), £(£) can be 
taken to be 1/3, since the amplitude of the oscillations 
remains unchanged for many oscillations. For the numer- 
ical values we have considered, we would get the produc- 
tion of helicity 1/2 gravitinos to be roughly 4 orders of 
magnitude larger than that of the helicity 3/2 states, if 
we naively assumed that helicity 3/2 saturated the Fermi- 
level (|/3| 2 = 1, which is an over estimated production). 
It is evident from the numerical example, that the pro- 
duction is suppressed by at least 8 orders of magnitude 
compared to that of helicity 1/2, see Fig. (|) and Fig. §). 



VI. IMPLICATIONS ON T RE h 

Through parametric resonance, other particles can also 
be created, and especially if their couplings are not sup- 
pressed by the Planck mass, they are perhaps produced 
more abundantly than gravitinos. Assuming that the 



end of reheating gives rise to a thermal bath with a fi- 
nal temperature T re h, it is possible to estimate the ratio 
n/s, where s is the entropy density, and n represents the 
number density of gravitinos after the end of reheating. 
Both n and s scales like a -3 , and we can also assume 
that the background scalar field density behaves on aver- 
age like matter, p^ oc a~ 3 , during the oscillatory period. 
It is possible then to estimate the final ratio by noticing 
that pi ps A 2 7Vg and n(t{) 
denotes the initial values 



A 3 -/Vq, where the subscript 



— rch ~ A7 - J- rch • 

s p t N 



(60) 



Here the left hand side represents the final abundance of 
gravitinos during nucleosynthesis. Gravitinos are weakly 
coupled to gauge bosons and its gaugino partners and 
their lifetime, T^ ccay w Mp/m^ 2 , is very long. For a TeV 

mass gravitino it could be around 10 4 — 10 5 seconds, and 
it would pose a genuine threat to nucleosynthesis. How- 
ever, this statement is strictly correct only for the helic- 
ity 3/2 component, since they can decay to gauge bosons 
and its gaugino partner through a dimension 5 operator. 
At high energies the interaction channels arc governed 
by 3/2 component rather than 1/2 component graviti- 
nos. In particular, helicity 3/2 may be produced with a 
mass close to a TeV range, so they decay very late and 
they are the ones which survive till late to cause prob- 
lems for nucleosynthesis. However, the same can not be 
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said with confident about the helicity 1/2 components. 
As we have seen, the violation of conformal invariance 
is not the same for both the helicities, and the helic- 
ity 1/2 gains an effective mass which is of the order of 
v$ 3> m 3 /2- Essentially, helicity 1/2 gravitinos are in an 
oscillatory scalar background with a frequency similar to 
their effective mass. There is no good reason to believe 
that the decay rate of helicity 1/2 to gauge bosons and 
to gauge fermions would mimic the decay rate similar to 
that in a flat background. So far, a detailed study is lack- 
ing in this area, but there is sufficient hint that the decay 
rate of helicity 1/2 is much smaller than the Hubble pa- 
rameter. We do not repeat the argument, rather we refer 
the reader to Ref. (lj|. The detailed calculation of the 
decay rate seems to be quite involved and we leave that 
for our future investigation. The important point to real- 
ize is that once the Universe reheats and thermalizes, the 
effective mass of the helicity 1/2 gravitinos becomes sim- 
ilar to that of the helicity 3/2, and as a result the decay 
rate would essentially be given by the usual decay rate 
in a flat background limit. Whatsoever be the detailed 



analysis, we must mention that while deriving Eq. (60) 



we have implicitly assumed that the initial abundance of 
gravitinos produced remain frozen until thermalization. 
Since the decay rate is smaller than the Hubble rate, the 
gravitinos produced during preheating will be able to sur- 
vive until the end of reheating. The final abundance will 
mainly depend on the model parameters such as A and 
Nq, see Eq. (|60|). Therefore, the constraint on the re- 
heat temperature derived from Eq. ( |60| ) is clearly model 
dependent. 

We also mention that we have not included the ef- 
fect of back-reaction coming from newly created bosons 
and fermions. In hybrid models, the production of the 
quanta associated with N and </> is very efficient and can 
take place just in the first few oscillations |10|, much be- 
fore the gravitinos have saturated the Fermi level. There- 
fore, back-reaction effects due to these quanta will quickly 
change the frequency and also the amplitude of the oscil- 
lating fields. We strongly suspect that especially in the 
hybrid scenario gravitino production will be affected due 
to such considerations and hence our current estimation 
of gravitino abundance will not hold anymore. However, 
even though the non-thermal production of gravitinos 
may stop after a while, there is a possibility to create 
them through scattering processes. This is beyond the 
scope of the present discussion and we shall hope to come 
back to these issues elsewhere [Bltl . 



VII. DISCUSSION AND CONCLUSIONS 

We have carried out the calculation for the gravitino 
production in multi-chiral field scenario, in particular in 
the context of hybrid model. As we have shown, it is 
possible to add a gauge-fixing term to the supergravity 
Lagrangian to get rid of the mixing between the goldstino 



and the gravitino field. Our method is a generalization of 
i?£ gauge studied in various contexts. We choose to work 
in a unitary gauge, in which the goldstino is completely 
removed from the physical spectrum. Our study empha- 
sizes major points in the non-perturbative production 
mechanism of gravitinos, analytically and numerically for 
the multi-chiral field models. We also give detailed analy- 
sis of the fermionic creation in general. We have observed 
that the fermionic creation in hybrid model is quite differ- 
ent from other chaotic inflationary models. The effective 
fermionic mass in hybrid models never vanishes, and as 
a result the particle production does not take place in 
first few oscillations, rather it builds up gradually. This 
makes it more interesting as far as the gravitino pro- 
duction is concerned. If we really want nucleosynthesis 
to be preserved in the context of supergravity inflation- 
ary models, we believe models based on hybrid inflation 
with low scales are probably going to be the only saviour. 
The reason is very simple; in other models there is no 
way we can argue the back-reaction due to the creation 
of other particles would stop creating gravitinos, but in 
hybrid models there is a scope where the back-reaction 
due to non-perturbative creation of bosons could affect 
the coherent oscillations of the inffaton and halt the par- 
ticle production completely. This gives us a new hope 
to understand the abundance of gravitinos during nu- 
cleosynthesis and we leave these important issues to be 
investigated in near future. 
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APPENDIX: ESTIMATION OF Kmak 

In this section we briefly discuss non-perturbative pro- 
duction of massive fermions. This discussion is general 
and self sufficient. We begin with the second order Dirac 
equation represented in terms of the two mode functions. 
It reads in conformal time 



(k 2 + (ma) 2 ± i(ma)') 



0. 



(Al) 



where u± are the mode functions, k is the momentum, 
prime denotes the derivative with respect to the confor- 
mal time, a denotes the scale factor, and m is an effective 
time-dependent mass. This mass can be written in terms 
of the physical time as 



m(t) 



m x + 9<i>{t) , 
<&(t)cos(m ( j > t) 



(A2) 
(A3) 
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where mx is the mass of the fermion we are interested 
in, 4>{t) is the scalar field oscillating with some initial am- 
plitude $(0), and m$ is its mass. The Yukawa coupling 
between the fermion and the scal ar is determined by g. 
Inspecting the mode equation Eq. ( |AlD , it is obvious that 
it mimics harmonic oscillator with a time-varying imag- 
inary frequency 



(ma) 2 ± i(ma)' . 



(A4) 



Particle creation occurs due to non-adiabatic evolution 
of the total frequency u>(t), that is, whenever 



duj(t) 
dt 



(A5) 



where we have expressed the total frequency in physical 
time. 

We shall explore here two possible scenarios for 
fermionic production, depending on whether or not the 
effective mass for the fermion vanishes at some point of 
the scalar field oscillations. To start our discussion, we 
consider first the situation when the amplitude of the 
scalar field initially satisfies |3>(0)| > mx/g, so that the 
effective mass goes through zero at some point during 
the oscillations. Since the particle production takes place 
within a time At <§; H^ 1 , we can safely neglect the ef- 
fect of expansion during the production of particles, espe- 
cially in a broad resonance regime^]. However, we have to 
keep in mind that the amplitude is gradually decreasing 
due to the effect of expansion. Once the amplitude drops 
below a critical value, the production of the fermions 
will soon stop and their number density will be frozen in 
time ra. The critical amplitude is given approximately 
by equating Eq.(A2) to zero 



4>* 



m x 
9 



(A6) 



Below this value, the effective mass remains always pos- 
itive. It is also important to mention that the produc- 
tion enhances near the regime when the effective mass 
vanishes, or, in other words, when there is a maximum 
violation of adiabaticity condition. With the above infor- 
mation it can be possible to estimate the typical momen- 
tum k required to violate the adiabatic condition when 
the amplitude of the field is close to This will oc- 
cur when cos(m^O ~ — mx / (g&*), or i» ~ ir/m^. Our 
condition Eq. (|A5|) implies 



^In the hybrid model scenario presented in this paper, it is 
evident from Fig.(Q) that the gravitino production is taking 
place in bursts within a short interval of time. However, there 
is another possibility, a narrow resonance regime, where the 
particle production takes place through out the evolution of 
the scalar field . 



2g 2 ^(t) 2 (m x + gcb(t)f + \g 2 '<j>(t) 2 > 

([k 2 + (m x +gmf) 2 +9 2 m^" 1 



(A7) 



Since, 4>(t*) = 0, and </>(£*) ~ mimx/g, we can estimate 



the left-hand side of Eq. ( A7) around 0(£*), and the final 
condition translates to a simpler form, 



1 



-m x m 6 > k 



3 

max • 



(A8) 




We have assumed uj ss A: max in the final derivation. This 
result confirms simila r res ult already obtained in Ref . [^| . 
We can express Eq. (A8) in terms of an effective q pa- 
rameter 



(A9) 



where q(t) — g 2 ^ 2 (t)/mi ) . In the above equation the 
value of q is evaluated at 4> ~ The q dependence in 
fcmax is quite different from the bosonic production. This 
is mainly due to the presence of the imaginary part of 
the frequency, which has a significant contribution to the 
violation of the adiabatic condition. At this point one 
may be able to estimate the maximum mass mx allowed 
by the violation of the adiabatic condition. With the help 
of the q parameter, it is possible to re-express m(i) as 



m(t) = mx + -^—^ — cosym^t) . 



(A10) 



where we have written explicitly the time-dependence of 
the amplitude. Hence, the maximum mass is achieved 
when = n/m^, and this gives 



7T 



(All) 



For reasonable values of the coupling constant g, it is 
possible to achieve very high values of the fermion mass 
mx ^> m^. This suggests that such production of 
supermassive fermions is indeed non-thermal and non- 
perturbative in nature. It is also important to notice 
that for values of q which are of the order of tens, the 
maximum fermionic mass obtained is of the order of the 
mass of the oscillating field. 

The above analysis suggests that creating super- 
massive fermions is possible because the effective mass 
m(t) vanishes around 0», and this is where the adiabatic 
condition is violated maximally. This happens quite nat- 
urally in quadratic inflationary potentials because the ini- 
tial amplitude of the oscillations for the inflaton is large 
enough to pass through the point where the effective mass 
m(t) vanishes A similar situation arises in the case of 
the helicity 3/2 states for the gravitino, and we find that 
its mass term Eq. ( |5l|) vanishes in each and every os- 
cillation of the N field, defined earlier. This means that 
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the adiabatic condition is violated maximally at those 
points. However, in the hybrid case the amplitude of the 
oscillations and the mass of the helicity 3/2 states are 
exactly the same mx = g$(i). This suggests that we 
can not create very massive 3/2 states, although we can 
create them with mass mx <C m§. Therefore, the hybrid 
situation is slightly different from the chaotic inflationary 
scenario with quadratic potential. 

Next, we study a scenario where the effective mass of 
the fermion never vanishes at any point, that is g$(0) < 
mx- This is the situation which arises in the hybrid 
inflationary scenario for the helicity 1/2 gravitino states. 
In this case the adiabatic condition is violated maximally 
near the point where m(t) is minimal but non zero ( 
then | <?</>(£) | is maximal), and this again happens when 
4>{t) ~ 0. Using Eq. (|A5| ), we can then easily estimate 
the upper limit on \<j){t)\ 



2/3 



V2 J 



- (m x - |#(i)l) 2 > k 2 



where we have replaced <j>(t) = —m 2 ^4)(t). Notice that we 
have explicitly taken the negative sign for <f>(t), in order 
to extremise the left-hand side of the above equation. 
We would like to see the range of the violation of the 
adiabatic condition for small k, 



mx 
3 



< 



\9<P(t)\ 



< 



m x 



(A13) 



Hence, in general the production could continue as 
long as the field amplitude follows the above condition 
Eq. Q. _ 

The maximal range of momenta for which the fermions 
are produced is obtained when |<7</>(i)| ~ mx, and is given 
by 



m\m x 



1/3 



V2 



(A14) 



and for mx ~ m^, the above expression reduces to 
fcmax < m x ■ This result could have been easily derived 



from Eq. (A8) by taking the masses to be almost equal. 

In the hybrid model, which we have considered, the 
amplitude of the oscillations is small, never exceeding 
mx/g- Hence we do not expect to produce super-massive 
fermions. The above result can be directly applied to 
the production mechanism of helicity 1/2 gravitinos. We 
have mentioned in the main section that the effective 
time- varying mass for the helicity 1/2 states does not 
vanish. The effective mass of the 1/2 states is almost 
equal to that of the oscillating frequency, which corre- 
sponds to having mx ~ m<h in our present discussion. In 
addition the adiabatic condition is broken for a narrow 
range of field values, see Eq. flA13| ). This means that the 
occupation number for the helicity 1/2 builds up in each 
and every oscillation and the gravitinos are produced in 
bursts. 
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